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We address quantitatively the relationship between the nonlinearity of a mechanical resonator and the non¬ 
classicality of its ground state. In particular, we analyze the nonclassical properties of the nonlinear Duffing 
oscillator (being driven or not) as a paradigmatic example of a nonlinear nanomechanical resonator. We first 
discuss how to quantify the nonlinearity of this system and then show that the nonclassicality of the ground 
state, as measured by the volume occupied by the negative part of the Wigner function, monotonically increases 
with the nonlinearity in all the working regimes addressed in our study. Our results show quantitatively that 
nonlinearity is a resource to create nonclassical states in mechanical systems. 

PACS numbers: 03.65.Ta, 03.65.Yz, 05.45.-a, 85.85.+j 


I. INTRODUCTION 

Mechanical systems are emerging as very well suited can¬ 
didates for the study of quantum behavior at the mesoscopic 
scale. Such a possibility is currently being explored in par¬ 
ticular in the domain of quantum opto-/electro-mechanics U), 
where ground-breaking experimental demonstration of quan¬ 
tum control of massive systems operating under explicitly ad¬ 
verse conditions have been recently made i). Yet, despite 
the substantial number of studies addressing the features of 
mechanical systems operating at the quantum level, only par¬ 
tial attention has so far been given to the physics of nonlinear 
mechanical devices. Classically, it is known that many non¬ 
linear systems exhibit very complex behaviors of potentially 
interesting features lO, and it is currently believed that such 
features might be potentially useful in many areas of investi¬ 
gation, even beyond physics. 

In the context of mesoscopic quantum behaviors, Katz et al. 
have studied the quantum-to-classical transition in the state of 
nonlinear nanoelectromechanical systems (NEMS) 0. Their 
analysis considered both an isolated resonator and one open 
to the effects of an environment, focusing on quantum sig¬ 
natures and on their disappearance toward classicality, as the 
operating temperature of the oscillator was raised. 

In this paper, we study the link between the enforcement 
of nonlinearity in a quantum mechanical oscillator and the 
manifestation of evidently nonclassical features in its state. 
Wondering about such a connection is indeed significant: the 
expectation values taken by observables of linear systems fol¬ 
low the corresponding classical equations of motion, and a lin¬ 
ear dynamics do not let phase-space non-classicality emerge 
(even at low temperatures). It is thus important, and relevant 
for applications, to understand which is the interplay between 
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the nonlinear character of the evolution of a given bosonic sys¬ 
tem and the strength of the nonclassical features that we are 
able to correspondingly enforce. 

From a quantitative viewpoint, here we will make use of the 
negativity of the Wigner function as a phase-space indicator of 
nonclassicality 0. This is an established notion of nonclassi¬ 
cality, with a close relationship with the non-local properties 
of the quantum state Ham. On the other hand, we will quan¬ 
tify the degree of nonlinearity of a given system by making 
use of the measures put forward by some of us in Ref. a. 
In order to complement the analysis presented in Ref. 0, we 
focus on the case of Duffing-type nonlinearity which, besides 
being technologically relevant as inherent in some forms of 
NEMS 0, has been the focus of a few studies on the quan¬ 
tum effects that it entails in terms of non-classicality noi and 
entanglement Ha. We show that a direct link exists between 
nonlinearity and the nonclassical character of the state of the 
oscillator. By focusing explicitly on the ground state of the 
system, we show that the phase-space nonclassicality of such 
state depends almost linearly on the degree of nonlinearity of 
the oscillator, thus suggesting a potential role of the latter fea¬ 
ture as a resource for the achievement of strong quantumness. 
Our conclusions are valid for a driven and an undriven Duff¬ 
ing oscillator, thus covering a vast range of physically relevant 
situations. We believe that our analysis and results embody a 
first interesting step towards the establishment of a rigorous 
link between such fundamental features in the dynamics of an 
oscillator. 

The remainder of the paper is structured as follows. In 
Sec. we discuss the specific example of nonlinear oscilla¬ 
tor, a Duffing oscillator, considered in this work. We focus 
on the undriven configuration of such oscillator and introduce 
the relevant measures of nonlinearity and nonclassicality that 
will be used throughout the paper. We show that a direct cor¬ 
respondence between degree of nonlinearity and nonclassical¬ 
ity can be established, pointing at the relevance of former for 
the enforcement of the latter in the ground state of an oscilla¬ 
tor. Sec.|n^deals with the case of a driven Duffing oscillator, 
and reports an analysis similar to the one presented in Sec.[n| 
We show that the relation between measures of nonlinearity 
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and non classicality is maintained in a dynamical situation as 
well, thus highlighting the fundamental nature of the relation¬ 
ship that we find, which appears to be um'elated to the details 


of the working conditions of the oscillator. Finally, Sec. IV 
closes the paper with some concluding remarks. 


with Af = 1/^1 + The corresponding probability den¬ 
sities are plotted on the right panel of Fig. [T] for increasing 
values of e, showing how the nonlinear term in the potential 
tends to localize the wave function of the oscillator around its 
equilibrium position. 


11. NONLINEARITY OF A DUFFING OSCILLATOR: 
UNDRIVEN CASE 


We consider a Duffing oscillator, which is described by the 
Hamiltonian m 

Hsys = - xj^cosojt, ( 1 ) 

where x and p are the dimensionless position- and 
momentum-like operators of the oscillator (such that [x,p] = 
i), £ is the anharmonicity parameter, is the amplitude of a 
possible force that drives the oscillator at frequency w. With¬ 
out loss of generality, throughout the manuscript we will con¬ 
sider the case of a stiffening nonlinearity with £ > 0. Such 
model is appropriate to describe the energy of a small-size 
doubly clamped mechanical resonator such as a carbon nan¬ 
otube, or a nanowire 0. The onset of nonlinear effects in 
such systems decreases with decreasing diameter of the de¬ 
vice, making either weak driving forces or thermomechanical 
noise sufficient to drive the motion away from the linear ap¬ 
proximation. In this Section we will focus on the undriven 
case, thereby setting 7^ = 0, deferring the treatment of a 
driven one to the next Section. 

It is instructive to gather an understanding of the form of the 
nonlinear potential energy to which the oscillator is subjected. 
In the left panel of Fig.fuwe show the function V(x) = £x^(4 
for different choices of the anharmonicity parameter, which 
shows that an increasing value of e results in more pronounced 
nonlinear effects at smaller displacements from the equilib¬ 
rium position a: = 0 of the oscillator. The effects of the quartic 
potential on the wave functions of the oscillator can be eval¬ 
uated by using time-independent perturbation theory. In such 
context, we use the notation Hsys = -f V(x) with 

= l{x^+f), V{x) = -^£x\ (2) 

We thus evaluate the first-order corrections to the eigenstates 
{|n)}ofijW as 


IV'n) 


k^n 


(3) 


We will mostly focus on the ground state (GS) of the nonlinear 
oscillator, for which we aim at finding an approximate form. 
While a fully numerical approach could be used to gather the 
full form of the GS of the oscillator, the range of values that e 
can take experimentally well justifies a perturbative approach. 
The perturbing Hamiltonian only couples |0) to |2) and |4), so 
that we have the normalised approximation of the GS 


IV'o) = A/" 




(4) 



FIG. 1. (Color online) The Duffing potential. In the left panel we 
show V(x) for e = 0.2 (red), e = 0.5 (green), e = 0.8 (blue). In 
the right panel we show the corresponding GS probability densities. 


In order to determine the range of values of £ within which 
the form of |'0o) given in Eq. 0 holds, we have calcu¬ 
lated numerically the GS lipnum) of the Hamiltonian Hsys 
in Eq. 0 using a truncated Hilbert space consisting of the 
first 51 number states and evaluated the state fidelity F{£) = 
|('0o|'0num)P. whose behavior against the nonlinearity pa¬ 
rameter £ is shown in Eig. Fidelity remains above 95% 
for £ G [0,0.8]. All the results reported in the remainder of 
this paper have been gathered using this range of values. As 
we show now, both the information on the modified potential 
energy and the GS wave function are key for the analysis at 
the core of this paper. 


F(£) 
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FIG. 2. (Color online) State fidelity F’(e) between the ground state 
of the undriven Duffing oscillator achieved through a perturbative 
approach |i/)o) and that estimated numerically by diagonalising the 
Hamiltonian model in Eq. 0 using the first 51 number states. 

We now pass to the introduction of the quantitative tools 
that we plan to adopt in order to gather insight into the rela¬ 
tion between nonlinearity and the onset of non classicality in 
the Duffing oscillator at hand. We would like to stress that the 
figures of merit that we introduce here go beyond the mere 
quantification of nonlinearity as given by the Hamiltonian pa¬ 
rameters (e.g., as given by e itself). This has the twofold ad¬ 
vantage of allowing to encompass situations in which more 
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than one Hamiltonian parameter is considered (see next Sec¬ 
tion) and of removing the dependence on the detailed form of 
the non-linear potential. 

We first consider a measure of nonlinearity based on the 
features of the GS ii and built by comparing ji/jg) and its 
unperturbed counterpart |0). Quantitatively, we determine the 
distance between the two GSs using the Bures measure: Given 
a perturbing nonlinear potential V(x), we define the nonlin¬ 
earity measure pB as the suitably normalized Bures dis¬ 
tance Dg between the GS of the oscillator under considera¬ 
tion and that of the corresponding harmonic one. In our case, 
we have 


Vb [V] = ^Dg Uo), |0)] = v'l - 1(0 |^o)|, (5) 

where we have used the fact that the two states under scrutiny 
are pure. As it is apparent from its very definition, this quan¬ 
tifier depends crucially upon the choice of a corresponding 
reference harmonic potential. Such a dependence can be over¬ 
come by considering a second way of quantifying nonlinear¬ 
ity: Given a potential V{x) with associated GS |0)y, we de¬ 
fine the measure of nonlinearity ?7 „g [V] 

?7«g[V] = ^„g[|0)(0|v], (6) 

where 5^0 [ 0 ] is the degree of non-Gaussianity introduced in 
Ref lfT3lfT^. This definition is intuitive: as commented ear¬ 
lier, a nonlinear potential would induce deviations from Gaus- 
sianity, which can in turn be used to quantify the strength of 
the nonlinear process itself. The degree of non-Gaussianity is 
built on the quantum relative entropy of the state and a refer¬ 
ence Gaussian state. As the GSs under scrutiny are always 
pure, we have 


me [V] = S[t] = h (^v^det[cr^]^ (7) 

with T the reference Gaussian state, S'[r] its Von Neumann 
entropy, (Tr its covariance matrix, and h{x) = (x -f |) ln(x -f 
i) —(x — ln(x— |). The crucial point here is that the defini¬ 
tion of Png requires the determination of a reference Gaussian 
state for the GS of V (x) rather than a reference harmonic po¬ 
tential for V (x) itself. Both measures are zero for a harmonic 
potential, whereas they may lead to different definitions of 
maximally nonlinear processes. 

While 0 < pb[V] < 1, the upper bound being reached 
if and only if the external potential V (x) gives rise to a GS 
orthogonal to that of the corresponding harmonic case, png is 
unbound from above, which complicates the quantitative com¬ 
parison between the two figures of merit. A suitable rescaling 
can be obtained at any fixed value of energy upon normalizing 
Png to the degree of non-Gaussianity of the states that achieve 
maximal non-Gaussian character at that value of the energy. 
This class includes number states and some specific superpo¬ 
sition of them (see m for details). The maximum of this 
rescaled quantity is thus achieved for a potential having a GS 
equal to a number state (n 0) of the harmonic oscillator or 
to some specific superpositions of them. 



FIG. 3. (Color online) Left: the nonlinearity measures ps [V] (red 
line) and pjvc [V] (blue line) for the undriven oscillator as a function 
of £. The inset is a parametric plot of ps as a function of pwe, show¬ 
ing that the two measures are monotone functions of each other. The 
right panel shows the Wigner function of the GS for e = 0.1 


One would expect that the nonlinearity increases with e. In¬ 
deed, as shown in Fig. this intuitive behaviour is captured 
by both measures, png [V] and pe [V], which grow continu¬ 
ously and smoothly with e. The two measure are also linked 
by a monotonic relationship. This is demonstrated in the inset 
of Fig.l^ where we show a parametric plot of p^ [V] against 
Vng [VjTthe curvilinear abscissa in such plot being embodied 
by e. Numerically, for e G [0, 0.8], the relation between the 
two measures of nonlinearity is very well approximated by the 
function 


Pb[V] =a + 6y/p7VG[V] (8) 

with a ~ 2.7 x 10“^ and b ~ 0.34. 

A broadly used indicator of nonclassicality in the state of an 
oscillator is provided by the volume occupied by its associated 
Wigner function in the negative region of the phase space 151 . 
The Wigner function of the state p of a single oscillator system 
is defined as 

Wpia) = ^j ( 9 ) 


where Xp(C) = Tr[pe^“^“^*“] is the Weyl characteristic func¬ 
tion and a G C. Unlike a true probability distribution, the 
Wigner function can take on negative values EMnl, which 
is a striking signature of nonclassicality. In the right panel of 
Fig. [^we show the Wigner function of the GS of the Duffing 
oscillator, showing the presence of regions of negativity that 
signal the nonclassical nature of the state of the system. In 
order to quantify such nonclassicality we use the measure 




Pp 

1 + 7p 


( 10 ) 


with p(p) = \Wp {a)\d'^a — 1 the negative volume of 

the Wigner function. The quantity p(p), which is per se suf¬ 
ficient to characterize phase-space nonclassicality, has been 
employed to study the quantum-to-classical transition in both 
linear and nonlinear oscillators iiiisi, as well as to charac¬ 
terize the performance of conditional schemes for the prepara¬ 
tion of nonclassical states of massive oscillators |l2l[T3- Here 
we consider its rescaled version according to Eq. (10 1 , which 
provides a number i/(p) G [0,1] that is thus amenable to a 
quantitative comparison with the proposed measures of non¬ 
linearity. 
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Making use of Eq. (j^ for the Wigner function of the GS 
of the Duffing resonator, we can evaluate the measure of non- 
classicality through a numerical integration. In Fig.j^we show 
the plot of the NG-based nonlinearity measure against the 
normalized measure of nonclassicality A numerical 

nonlinear fit gives the functional relation 

= 0.002 + 0.207yi^(V'o) + 2.731:/(V'o), (H) 

showing that after an initial trait, the link between the two 
quantities becomes linear, ensuring the proportionality of the 
two figures of merit under scrutiny. 

Vng 




FIG. 5. (Color online) In the left panel we show V‘^{x, t) against 
time and position for = 0.015, e = 0.1 (red, opaque surface), and 
e = 0.5 (blue, transparent surface). In the right panel we show the 
corresponding GS probability densities, \%l)o{x,t — l)pforthe same 
set of parameters. 


when the perturbative potential is off, as 

IV'(O)) = ^C„(0)|n) (c„(f) e C) (12) 

n 

and aim at finding a perturbative expansion Cn ^ (f > 0) S C at 
order q G Zin the perturbation, so that the state of the system 
at the corresponding order becomes 

(13) 


FIG. 4. (Color online) We plot the the nonlinearity measure [V] 
against the normalized measure of nonclassicality based on 

the negative parts of the Wigner-function for the GS of the undriven 
Duffing resonator. We show that the two measures are monotone 
functions of each other. The squares show the values of rjNG at set 
nonclassicality of the ground state of the undriven Duffing oscillator. 
The solid line represents the fit given by Eq. fn). After an initial 
trait where ting grows as the relation between the two 

figures of merit becomes approximately linear, showing the direct 
connection between nonlinearity and nonclassical character of the 
state of the oscillator. 


III. NONLINEARITY OF A DUFFING OSCILLATOR: 
DRIVEN CASE 


We shall now pass to the analysis of a driven Duffing 
resonator, which is an example of nonlinear resonator often 
encountered in relevant experimental situations 
The dimensionless Hamiltonian of the system is thus Eq. 0 
with the explicit inclusion of the time-dependent driving term 
—xJ^ cos ujt. In what follows, we choose a working point well 
within the region of bistability of the oscillator (which is en¬ 
sured for G [0.015,0.06] and w G [1.016,1.02]). The form 
of the driving potential V‘^{x,t) = — xTcos ujt is il¬ 

lustrated in Fig. 1^ showing the deformation induced by the 
nonlinear and time-dependent part of the perturbation. In or¬ 
der to evaluate the form of the GS associated with the full 
model, we shall resort to the use of time-dependent perturba¬ 
tion theory. We decompose the state of the system at f = 0, 


where is the eigenvalue of the unperturbed Hamiltonian 
corresponding to the eigenstates jn). The first-order correc¬ 
tion in both E and e, which will be the highest order of the 
perturbative expansion that we will consider here, is given 
by El 


c^uHt) = -^ ( 14 ) 

^0 

with = 6ni, = {n\V‘^{x,t)\l) and ^ 

The explicit calculation of [^/’(f)) up to the stated order 
of approximation and for [■(/'(O)) = [0) leads us to the GS 
wavefunction 




+E 

e 

^16 


(1 — e®* [cos ojt — iuj sin cot ]) _ 3 , 

V2(a;2 - 1) ® 

372(6-5®*/^ - e-®*/2)V'2(cc) 




/'g 2^* — e 2^') 

+ ^^- h^ix) 


being A4 the normalisation factor and tpn{x) = {x\n) the 
wave function of state \n). A plot of the spatial distribution 
\ipo{x, f)|^ at a set time t is shown in the right panel Fig.|^ 
Looking at the potential, one would expect that the nonlin¬ 
earity increases with e at any fixed values of the given param¬ 
eters. Indeed, this intuitive behavior is captured by both mea¬ 
sures, 77 „g [V] and [V], as they grow continuously. The two 
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FIG. 6. (Color online) Left: the nonlinearity measures rjB [V] (red 
line) and [V] (blue line) for the driven oscillator as a function of 
e and for the choice of parameters T — 0.015, t = 1 and to = 1.018 
which ensures that the resonator is in the bistability region. The inset 
is a parametric plot of r/B as a function of rj^G, showing that the 
two measures are monotone functions of each other. The right panel 
shows the Wigner function for the same set of parameters and e = 
0 . 1 . 


measures are monotonic functions of each other, as illustrated 
in the left panel of Fig. In the right panel we also show the 
Wigner function for t = 1, which shows negative regions and 
thus the signature of nonclassicality. When assessed against 
the chosen indicator of nonclassicality, ijnc is again found to 
be in direct cori'espondence with the volume occupied by the 
Wigner function corresponding to |'0o(f)) [cf. Fig.|^. 


Vng 



IV. CONCLUSIONS 

In this paper, upon using purpose-tailored quantitative fig¬ 
ures of merit, we have addressed in some details the relation 
between nonlinearity and nonclassicality in a class of nonlin¬ 
ear oscillators that is relevant in many contexts, including ex¬ 
perimental solid state physics. By approximating the form of 
the ground state of the oscillator through a perturbative ap¬ 
proach (either stationary or time-dependent) we have been 
able to demonstrate that the negativity of Wigner function, 
which is a well-acquired measure of nonclassicality in contin¬ 
uous variable systems, is in monotonic relation with recently 
proposed measures of nonlinearity, therefore reinforcing the 
idea of nonlinearity as a catalyst of quantumness. Although 
our conclusions have been gathered by addressing the specific 
example of a Duffing oscillator, our work paves the way to in¬ 
teresting extensions, primarily concerned with the application 
of our tools to other forms of nonlinearities M- 

A second interesting direction of investigation would deal 
with the inclusion of environmental effects, and with the pos¬ 
sibility to shield the degree of nonclassicality enforced in the 
state of a quantum oscillator through a suitable degree of non¬ 
linearity. This might entail an interesting way of protecting 
quantumness, stemming from the direct, non-demanding con¬ 
trol of the Hamiltonian of the oscillator. In fact, while the har¬ 
monic assumption is an approximation valid within many con¬ 
texts (from nano-mechanical oscillators to ultracold atomic 
systems in external potential), switching to explicitly non¬ 
harmonic situations is typically straightforward by the means 
of a strong driving. This is generally more economic than 
time-gated external pulses (required in dome of the techniques 
devised so far for the protection of quantum features) or the 
control of the properties of the environment, which is typically 
of not easy access. Work along these lines is in progress and 
results will be presented elsewhere. 
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